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Abstract 

I The colored objects - quarks and gluons - being confined in a small volume ^ ~ i?Q, 

^ ' Rq ~ 0.5fm inside the QCD bound state get there not small masses mqq ~ IGeV, 

fSj ' rrig ~ 0.5GeV. This drastically simplifies the QCD dynamics, as now the probabilities 

e.g. of production of one extra massive valent gluon or extra qq pair, turned to be 
small due to a large gap between corresponding energy levels. The ordinary quantum 
mechanical perturbation theory calculations made in the paper shows that corresponding 
. dimensionless transition amplitudes Amh = — (of'^ (o) (from meson qq to the hybrid 

qqg) have a value of 20-25% and the mixing of hybrid with glueball (gg) is even smaller 
(of the order of 10-18%). Taking into account all such mixings with all highly lying 
Fock states is equivalent to the exact solution of QCD bound state problem. This can be 

■ done really summing the small perturbation theory results for these states contributions. 

o 
o 

O ■ 1- INTRODUCTION/MOTIVATION 

^! • Massless quarks and gluons acquire inside of hadrons not small constituent masses 

^ : mg ~ — ~ 0.35GeV, nig ~ 0.45GeV 

■ due to gluon string confinement in a small volume Vq ~ Rq, Rq — 0.5 — 0.6fm. Also the 
. gluon strings themselves are massive having about the same mass M string ~ 0.4 — O.SGeV, so 

5^ ! mg ~ Mstring at String tension of about 0.18-0.15GeV^. 

• Therefore each extra gluon, or a qq pair of quarks lead to enhancement of hadron mass 
to about 0.8 — l.OGeV. This was confirmed by lattice calculations. 

• These estimations can be illustrated comparing e.g. the p -meson mass nip ~ 
O.SGeV (p = (qq)^^^) with the mass of vector hydrid h = {qq,g)'^^^) for which nih — 1.8 
GeV. Another example represents glueballs, G = {gg), consisting from two valent (con- 
stituent) gluons bound by a pair of gluon strings and observed apparently with the mass 
mo — 1.5 — l.SGeV. These recent data agree well with the estimations given above and with 
lattice calculation data. 

• This situation can naturally be understood dividing the gluon field in two parts (as 
t'Hooft and Simonov [1] have suggested): A^(x) = B^{x) + d^(a;) , where B^{x) = taB'^{x) 
leads to gluon string formation and the part 

a,{x) = t^a-[x) = V'^[c;(^)e+^^^-^ + c^+(fc)e-^^-1 (1) 



is responsible for the valence (constituent) gluon production and absorbtion. 



The QCD interaction energy has the form 

Vqcd = VQcoia^,) + Vqcd{B^,) + V'^^^^b^- 

where term 

^qcdI^m) =3s j [{q{x)li,tc,q{x))al{x) - ]^fc,p^ala^^{dxal - d^al) - 0{al)]d^x (2) 
changes the numbers n, of valence gluons and of qq pairs respectively; the term VQ^^{Bfj,) 

3 

leads to QCD strings production and to colour confinement, while K' ^ = Cfc(o./u)fc(i?/u)^~'^ 
is not important below. 
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2. THE ZERO ORDER HAMILTONIAN AND BOUND STATES IN THE 
SEPARATED QCD SECTORS. 



Considering v'q(jj^ as perturbation let us introduce the zero order QCD Hamiltonian 



^QCD- 



Hqcd = Hqcd + ^qcd(«m) (3) 
with the complete and orthogonal system of its eigenf unctions 

^QCD^N — \^) 

of separate Fock sectors with a fixed numbers n, of gluons and of qq pairs (at fgc^) = 
their numbers n, N remain fixed!) Both the QCD Hamiltonian Hqcd (3) and its zero order 
formQ Hq(jj^ were introduced and intensively used in papers by Simonov and collaborators 
[1,2, 3]. Their results are largerly used below. Let us stress that having the values of all mixing 
rates \xf\'^ i.e. of all Cf^ = CjXf ,one directly recovers using Eq.(5) below the exact QCD 

solution ^! j of bound state problem. It is explained farther that A^'^'' can be found in the form 
of usual quantum mechanical theory rapidly convergent series considering QCD interaction 
'^'qcd of constituent quarks and gluons in Eg. (3) as small perturbation. 

• The exact solution ^ j of QCD bound states problem can be written in general as super- 
position of the zero order eigenf unctions ^j^'' of Hamiltonian ifg^^: 

oo 

= E CN,n^N\ With \CiN? = (^;^.) = 1 (5) 

N,n=0 n,N 

as (<l>i;^<l>y*) = 6n,n'SN,N'- Here j- can be the meson j = M (qq), the glueball j = G = (gg), 
or different hybrids: j = H = (qqg), or j = H' = {qq2g) etc. 



It is usefuU to put the superpositions (5) for j = M, G, H in the form of Fock columns. 



which includes only the string interaction part Vq(^jj{Bij^) of the total QCD interaction. 



Denoting N,n — k — (M, G, H), C' ' — CjXj ' one has respectively: 
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(6) 

where Cm = (l - l^Z? - {X^m^I^Y^^ ^ I ^ Co C,, and a\\ \f' < < 1 as the mixings 
with higher Fock states A 
below. 



(k) 
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(0) .,(0) 



■C 1 are proved to be small (A^'^'' ~ 10 — 20%) - see 



3. THE QUANTUM MECHAMICAL (QM) PERTURBATION APPROACH 
(PA). 

• Now we are ready to apply the usual QM perturbation theory (PT) which gives the well 
known results: 



+ E 



n 



TTT + ■■■ 



^ r 



^n[l - \ E 



where 



Vkn = J K^QCD^kdr 



0T + 



(7) 
(8) 



and operator Vqcqo transists $m so that 14'^^ 7^ only for k ^ n 

• Let us put the zero order wave functions for j^'^ — O"*""*" and 1 mesons in the form: 



(9) 



separating vertices R^'^ of meson decays into quarks and gluons (shown in Fig. 1,2) from the 
radial parts 954^(^12)5 V^i?(^i25 "^13)1 VGifgg) of the corresponding wave functions; here a,b = 
1, 2, 8, A = 1, 2, 3 are gluon and quarks colour states indices wile /x, = 1, 2, 3, 4 are gluon 
Lorents indices and X, (3 = 1, 2, 3, 4 are Dirac quarks indices. 
The explicit forms of vertices for 0+ and 1~ cases are: 



( (f)0+\\a _ 
\^M)al3 - 

/ dO+ \Ao-,a _ 1 +\a^a(^, \ 

)a/3,ij. ~ 4V2 « '/^^"/^ 

/r>0+\ab _ _J_Xab„ 



(^^")aAv = i^fe;i(7.7M)a/3 



(10) 



they are normalized to 1 : {R\R-'2) = 1 for all A = M,H,G, where a = 1,2, ...,8 and 
ia = are Gell-Mann 3x3 matrices. 

• With these zero order wave functions the QCD perturbation operator v'qq]^ (see the 
first term in Eq.(2)) can be used in calculation of Vkn directly in r— space in the following 
"stripped" form: 

{gs/V2^c)ilMu)[S^'\n3) + S^'\f23)] (10a) 



as quarks g, q and gluons g individual wave function can be absorbed (in Eq. (8) for Vkn) 
into $fc and while (5*^^-' terms prove to appear in {vQQD)stript as e.g. transition M ^ H is 
reached by production of a gluon by quark or anti quark at some point f = and also by 
absorbtion and production of quark essentially at the same point r = fi, i.e. at ri3 = (or 
at the point r = r2, i.e. at r23 = for anti quark). For H G transition the situation is 
simpler and here the action of fgc/) involves the appearing of 5^^\fi2): 

due to a locality of the QCD interaction; see the first term in Eq.(2) (acting on the hybrid 
state it absorbs q, q in some point x = Xi and produce gluon g just in the same point). 

This "stripped" version of PT follows simply from the momentum representation (like one 
used by Orsay group paper [4]) of all constituent gluons and quarks operators, e.g. for quarks 

q^{x) = ^^-^[«„(p,a)^,e^^^ + <(p,a)t;,V^^] 

and similarly for gluons by the above determined operator a^{x). 

Authors of two recent publication [6,7] (from North Carolina Univ.: E.Gubankova, C.-R. Ji 
and S.R.Cotanch) have calculated the glueball spectrum and obtained the value of gluon mass 
rUg ~ 0.5GeV inside of confinement region - close to the Simonov results [1] used here. They 
proposed to consider the colour Coulomb interaction of quark and gluons as the sours of mixing 
between different states in each QCD isolated sectors (with a given numbers Ug of gluons and 
N = Ngg of valent qq pairs). Let us remind that t'Hooft identity {A^{x) = B^{x) +a^{x), see 

above) allowed us to put the QCD interaction in the form Vqcd = ^gcni^iJ.) + ^qcd('^m) "I ' 

with Vq(jj~,{B^) responsible for the binding of all separate Fock state (with a given ra^, 
see Eq.4). 

However, a short disqussion made above has clearly shown that, while the Coulomb coulor 
interactions of constituents inside of each sectors can be important as authors of Refs.[6.7] 
notes (less important, however than the string type one), the real cause for mixings of bound 
states is the QCD interaction (2) v'qfjj^ia^) of valence gluons and quarks. It changes the 
numbers of them converting e.g. the mesons {Ug = 0, ngg = 1) into hybrid {ug = Ugg = 1) and 
in glueballs {ug = 2, Ngg = 0) and really mixing them. Just the rates of these mixings are 
calculated below. 

• The correction (7) to En and to \E'„ of PT determined above proves to be not too small-of 
the order of 10-20%, LcA^f ^ = < 0.10 - 0.20, for A,B = M, H, G. Therefore one can 

calculate them using approximate zero order wave functions. The explicit form of meson and 

^Physically it means that gluon strings stretched between valent quarks or gluons can have the transversal 
oscillations. After quantisation the quants of these oscillations will represent the valence gluon states. 

The potential WqcdI'^a) determines their interaction whith quarks and leads to production or absorption of 
these valence states of gluons and of qq pairs. 



A 



glueball wave functions can be determined from Schroedinger equation with the string type 
Unear rising potential Vstring — (^l^i — O.lSGeV^, of qq interaction or respectively, by the 
potential Vgg ~ (9/4)0-17^ for gluon-gluon interaction case. They can be approximated (with 
high accuracy) by a simple Gaussians: 

(11) 

gO 

whith the proper chosen values of Tq and r^Q which are 

rg = {f^)j^ = 0.725/m^ Vo ~ A.32(GeV)-^ 



r 



2 _ 
9 



respectively obtained at cr = 0.18GeV^,and at agg — jcr. 

• The hyperspherical formalism [4,1] allows one to obtain with the same accuracy the 
three-particle (gi, q2 and g) hybrid's wave function (7^12^3) = p~^/'^XH{p)U'k{^a) , where 
Uk{^a) = const for the hybrid ground state, and 

= + if = E '^r^k with /i = /ii + /i2 + /is = 2/ii + /is at /Xi = /X2 = /Xg, /is = l^g- 
i,k 

At ri3 = one finds fi^ = fi2 and = ^^^2mM^^^ ^3 '^h®^® /i + /ia = 2(/ii + //a) and the scale 
m can be choosen arbitrary. Putting /i • mo — /ii(/ii + /is) or nio — /ii ^i^^+2^\ one obtains at 
ris = 0: [p\i3=o = r^g = = r^. 




Note that at ri2 = ^ (^^^)i/2r23 (here rig = ras), i.e. p = Oor^ with 
ao = [2/(1 + /ii//is)^]'/' ^ 1.083 at /ii = /i, = 0.312GeV, /ig = 0.44GeV, m„ = /ii^g^ 
(in the glueball /ia = (/ig)G = 0.53GeV, while inside of hybrids [1] /ia = (/ig)i? — 0.44GeV). 

In the hyperspherical method the total potential [4,1] 

V^(p) = x„.p+^|±i> (12) 

(where L = | + /sT, K ^ 0,1, L = | at = and xq = ryo^^ = t^q ]^^'^^' , 

VMM -y/l+2|Ui//i3 

770 = -^5/2 = 1.357) has a minimum at p = po ~ l/(m(7)^/^ where VF'(p) = 0, i.e. at: 
xoCT = (see Fig.5.), resulting in: 



\pm \i\ + /is xqc^ 



1/3 



(13) 



Therefore for the lower hybrid state one obtains practi- 
cally the oscillation potential (see Fig.5): 



with .i = w'(p.) = m^, -0 = ^ (") 

Ail = A*2 = A*g, A*3 = "^s, A* = 2mg + 

This gives the normalized hybrid wave function in the 
following simple Gaussian form similar to Eq.(lO): 




951/(^3^23) 



\det\ 



-1/2 



.(!!^)l/4^-5/2e-^(p-P0)^ 



Fig. 5 



(15) 



where — {-^)\y + |(1 — ■£;) + ■■■] will be used for the effectively small Ap — p — Po Po, 



PO' 



\det\y-^\det[ % 'I |')|V2=f/W3\ 
i T] R ) \ y 



3/4 



while = = . 3, 

and the normalization means 

^/dVidVadVsl^ifp = /d=^CA|dei||0j^p = / \det\p^dpd^^\(^H\^ = 1 



4. MATRIX ELEMENTS OF MH,HG,MG; MIXINGS AND CORRECTIONS 
TO PARTICLE BOUND STATES ENERGIES (MASSES) 

• Matrix elements Vuh-, Vhg, Vmg are generically small - of order ~ (1/5) — (1/10) GeV; 
let us show this for the simplest 0^"^ mixings case: 

a) The meson-hybrid mixing amplitude is represented in the first PT approxima- 
tion by two graphs 3a,3b,yielding: 



as 



Vmh = 2-^ J ^^(r)(^^(7.e^4)^^^)(^^(r)dV = 

(^^"(7.e^4)0 = ^. 
Vmh = AMHJ{a,b) 



3 



Thus 
where 



A 



(16) 
(^•7\ 



while 



oo 



J(a, b) = 1/6^/2 / g-t2_„(,_,)2^2^ + -(1 - t/b) + ...], t=^ (U 

2 2ro 







and the decomposition of (fn in power of = ^ <C 1 (shown above after Eq.(15)) was 
used, so that (^^(ria = r, = 0) ~ f-5/2 ^-a(t-b)'^ written at \t — b\ 1/a < b as: 
^-5/2g-a(t-6)2 + 1(1 - t/b) + f (1 - + ...] with dimensionless a = ^^r^, 6 = 



with po; = (y/^i + ^o)e// determined in Eqs.(13),(14) through quark and gluon masses 
/^i, /i3, resulting in po — 3.86GeV~^, cuq — 1.02GeV. 

The simple numerical calculations yield: a ~ 2.79, b ~ 0.775, and Amh — 0.694GeV 
obtained at Qs = 2.24, ag ^ 0.4, while for the rapidly convergent integral (18) (where the 
integrand in J(a, b) is shown in Fig.6a) one obtains: J(a, b) — 0.290 in the first approximation. 
This results in: 

Vmh = ^Mi/ J(2.79, 0.775) ~ 0.694GeV • 0.290 = 0.201GeV, 

leading to about 22% for MH dimensionless mixing amplitude: Xmh = i pol'lo i — 0.22 for 
\E^-E%\ -O.SGeV". 

b) Quite similarly for hybrid -glueball mixing one has (see two graphs in Fig.4, 
each giving the half of Vhg value ) : 

2 V 8 Tg TT-^/^rg 1 + 2/X2/At3 <^o 

and 

J = J(a, b) = b-'l- J e-^'-'(^-''n^dt[l + ^(1 - ^) + f (1 - + t = ^aor (20) 







as here {R^H{luelk)R^G) = ^nd one had p = cor, cq = sjj^^ - 1-083 at ri2 = 0. 

Also: as ^ ~ 0.719 then: 

~ _ 2r2^^^^^j2 _ (^fa ~ 0.517a ~ 1.444 
6 = ^ = 6/(^) ~ 6/0.719 = 1.078 

Simple calculations of J = | J(a,6)|(see Fig 5.b) gives: J ~ 0.285, and thus: 

Vhg = AhgJ = 0.320GeV • 0.285 = O.OgiGeV 

with the first correction in Eq.(20) at Ahg = 0.320. That results to about 18% of H-G mixing 
amphtude (at AEhg = - E'^^l - 0.5GeV): 

Ai/G = 1:^,^^1-0.18 -18%, 



and to very small amplitude Xmg of MG mixing: 

Amg = Xmh^hg ^ 0.045 ~ 4.5% 

• The value of AE^ = Am„ = ^ jX"^po in Eq- (7) gives in CM system of any particle 

the correction to its mass m° due to the mixings with higher Fock states. These corrections 
to meson, to hybrid and to glueball masses prove to be very small: 

AruM ^ -AruH = -^^f ^ -50MeV, 

where the energy gups AE%jfj and AE'jjj^ were put for estimations to be equal to some fraction 
of IGeVXAE^QH ~ 0.5GeF, AEmh ^ 0.8 GeV) 

There exists a danger that highly exiting meson states with Ej^ ~ IGeV can enhance 
significantly these estimation for Auih- 

5. Summary. 

• It is argued that gluon field can be represented as the sum of continuous part B^{x) 
leading to gluon string creation and to confinement of valent (constituent) quarks and gluons 
inside hadrons in a small volume V ~ i?Q, Rq ~ 0.5 — 0.6fm and the part a^^x) which 
produce and absorb these valent quarks and gluons. Due to the confinement these valent 
colour particles acquire not small mass ~ l/i?o ~ 0.3 — 0.4GeV. Also gluon string itself has 
about the same mass. 

• Therefore each extra gluon, or pair of quarks qq bound by gluon strings lead to enhance- 
ment of hadron mass to about of \AE\ ~ (0.8 — 1.0)GeV 

• In the quantum mechanical (QM) perturbation theory (PT) this large energy gap leads 
to rapidly convergent series as the PT matrix elements of QCD interaction (2) ^^^^(a^) = 
9s J{(llfitaq)a'^{x)d^x connecting these sectors (like meson qq, hybrid qqg, glueball gg with 
different number of valent gluons and qq paires) are small: Vmh — 0.2GeV, Vhg — O.lGeV, 
Vmg - 0.02GeV. 

These figures and those given above were obtained in the first approximation in expansion 
of integrand in Eq.(18) in powcs of ((1— It converges only asymtotically (up to n < 3—4) 
and the odd powers of n gives to J(a, b) the negative contribution. The next two terms in 
Eq.(18) up to n < 3 enhance slightly the value of MH and HG mixings changing nothing 
in princiole (as Vmh changes from ~ 0.20GeV only up to 0.26GeV while Vhg encreases from 
O.OQGeV forn < 1 up to 0.013GeV at n < 3 in Eqs.(18)-(20). 

• As a result the QCD problem of colourless bound states can have a simple exact solution 
in the form of convergent usual QM perturbation series which converges rapidly just due to 
the confinement effect. 

• The developed logic can have a few dangerous underwater stones: in the suggested 
"stript" relativistic approach the small dependence of functions ^M{f), ^G{f), ^if(^, ^3) on 
quark spin variables was not taken accurately into account. The more accurate approach is 
now in progress. 



n 
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